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We put constraints on primordial black holes (PBHs) as dark matter candidates by considering 
their capture by neutron stars (NSs). If a PBH is captured by a NS, the latter is accreted onto the 
PBH in a very short time leading to the star's destruction. Mere observations of NSs put limits on 
the abundances of PBHs. We are able to exclude PBHs as dark matter from observations of NSs in 
globular clusters for PBHs masses 10 18 g < ibbh < 5 x 10 23 g. 



I. INTRODUCTION 

The existence of the dark matter (DM) has been es- 
tablished so far only through its gravitational interaction. 
Consequently, little is known about the DM nature apart 
from the fact that it is non-baryonic, non-relativistic, 
weakly interacting and constitutes about 23% of the total 
energy budget of the Universe (for a recent review see, 

e.g., BE!). 

Various candidates for the DM have been considered 
in the literature. In the context of particle physics they 
are associated with a new stable particle beyond the 
Standard Model, a popular example being the so-called 
Weakly Interacting Massive Particles (WIMPs). How- 
ever, candidates that do not require new stable particles 
also exist and are still viable. An attractive candidate of 
this type is primordial black holes (PBHs) O H] . This is 
the possibility we consider in this paper. 

In the early universe, some primordial density fluctua- 
tions could have collapsed producing a certain amount 
of black holes. These PBHs possess properties that 
make them viable DM candidates: they are nonrela- 
tivistic and have a microscopic size of the order r ~ 
10 _8 cm (mBH/10 20 g), which makes them effectively col- 
lisionless. The initial mass function of PBHs depends on 
their production mechanism in the early universe and is, 
essentially, arbitrary. 

There exist a number of observational constraints on 
the fraction of PBHs in the total amount of DM. First, 
PBHs with masses tobh < 5 x 10 14 g evaporate due to 
Hawking radiation [5] in a time shorter than the age of 
the Universe and cannot survive until today. At slightly 
larger masses, even though the PBH lifetime is long 
enough, the Hawking evaporation still poses a problem: 
the PBHs emit 7-rays with energies around lOOMeV [6 
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in the amount that contradicts the data on the extra- 
galactic gamma-ray background. For instance, the Ener- 
getic Gamma Ray Experiment Telescope j7| has put an 
upper limit on the cosmological density r2pBH < 10~ 9 
for ttibh = 10 15 g [5]. From such observations, one can 
infer that PBHs with masses msH < 10 16 g cannot con- 
stitute more than 1% of the DM. More massive PBHs 
were constrained by EROS microlensing survey and the 
MACHO collaboration, which set an upper limit of 3% 
on the fraction of PBHs in the Galactic halo in the mass 
range 10 26 g < m B H < 10 30 g \W\. Theoretical pre- 
dictions on the abundance of PBHs have also been com- 
puted using microlensing of Kepler stars [TT1 IT2"] . At 
even larger masses 10 33 g < tobh < 10 40 g, the three- 
year Wilkinson Microwave Anisotropy Probe (WMAP3) 
data and the COBE Far Infrared Absolute Spectropho- 
tometer (FIRAS) data have been used to put limits 
on the abundance of PBHs [T3|. These constraints are 
shown in Fig. [TJ They leave open the window of masses 
(a few) x 10 16 g < m BH < 10 26 g. 

In order to put constraints on PBHs in the remaining 
allowed mass range, in Ref. |14| we have considered the 
capture of PBHs by a star during star formation process 
and their further inheritance by the star's compact rem- 
nant, the neutron star (NS) or the white dwarf (WD). 
The presence of even a single PBH of a corresponding 
mass inside the remnant (NS or WD) leads to a rapid 
destruction of the latter by the accretion of the star mat- 
ter onto the PBH [T5rfT9] . Thus, mere observations of 
NSs and WDs in a DM-rich environment, such as centers 
of globular clusters, impose constraints on the fraction of 
PBHs in the DM and exclude PBHs as the only DM can- 
didate in the range of masses 10 16 g < ttibh < 3 x 10 22 g. 
Still, a range of PBHs masses from 3 x 10 22 g to 10 26 g 
remains unconstrained. 

In this paper we derive the constraints that arise from 
the direct capture of PBHs by already existing NSs. The 
origin of the constraints is the same as in Ref. [H] : even a 
single PBH captured by a compact star rapidly destroys 
the latter, so the existing observations of the NSs and 
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WDs require that the probability of capture is much less 
than one. Assuming that, we put constraints on the PBH 
abundance. 

We find that the argument based on the capture of 
PBHs by the existing NSs allows one to extend the con- 
straints to larger PBH masses and exclude PBHs as com- 
prising 100% of the DM up to m B H < 5 x 10 23 g, leaving 
open only a small window of about two orders of mag- 
nitude. Also, the constraints on the fraction Opbh/^dm 
of PBHs in the total amount of DM at large PBH masses 
become tighter as compared to Ref. [TJ]. The final situ- 
ation is summarized in Fig. [T] 

The rest of this paper is organized as follows. In Sect. 
II we discuss the capture of PBHs by compact remnants 
due to dynamical friction. In Section III we derive the 
constraints on the fraction of PBHs in the DM from the 
observations of NSs. In Section IV we summarize the re- 
sults and present our conclusions. Throughout the paper, 
we use the units h = c = 1 



II. CAPTURE OF BLACK HOLES BY 
COMPACT STARS 

A. Energy Loss 

A PBH is captured if, during its passage through a 
star, it loses its initial energy and becomes gravitation- 
ally bound. From this moment every subsequent PBH 
orbit will again pass through the star, so that finally 
the PBH will lose all its energy and will remain inside 
the star all the time. Therefore, the criterion of capture 
of a PBH is i^ioss > with E\ osa being the en- 

ergy loss during the collision and v the PBH asymptotic 
velocity. Two mechanisms of energy loss are operating 
during the collision: deceleration of the PBH due to the 
accretion of star's material and the so-called dynamical 
friction [H] . In the relevant mass range of PBHs the 
former mechanism is less efficient compared to the latter. 
We, therefore, focus here on the energy loss due to the 
dynamical friction. 

As a PBH passes through the star, it transfers momen- 
tum and energy to the surrounding matter. The result, 
called the dynamical friction, is a net force that is oppo- 
site to the direction of motion of the PBH. As long as 
the PBH velocity v during the collision is larger than the 
velocity of the particles constituting the compact object 
(which is a good approximation for compact stars), one 
may take the dynamical friction force to be 
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where p is the density of the star matter. For the sake 
of the estimate, it can approximately be considered as 
constant, p — 3M /(4-7ri? 3 ), with M and R being the mass 
and the radius of the compact star. The factor ln(A) is 
the so-called Coulomb logarithm [201 HI] • 



The total BH energy loss can easily be computed from 
cq. . Taking the PBHs velocity during the collision to 
be v = v csc — yj2GM/R >• Vo, the energy loss during 
one collision is 



-^loss — 



R 



(2) 



As one can see, the energy loss depends on two param- 
eters: the PBHs mass and the radius of the compact 
object. Since E\ oss is inversely proportional to the ra- 
dius of the star, NSs induce a larger energy loss during 
one collision compared with WDs. Thus, we will only 
consider the case of NSs from now on. In this case 
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where we have substituted R = 10 km. 

A remark is in order. The core of a neutron star is 
comprised of the degenerate neutron gas, so the question 
arises whether eq. [T] is still applicabk£] Note, however, 
that by the time the falling PBH reaches the star it picks 
a relativistic velocity v ~ 0.6c and can transfer to neu- 
trons the momentum of up to ~ 1.8 GeV, which is by a 
factor of a few larger than the Fermi momentum of neu- 
trons in the center of the star, and much larger than the 
Fermi momentum away from the center. So, while some 
reduction as compared to the estimate ^ is expected, 
this reduction is not dramatic and is comparable in mag- 
nitude with other uncertainties involved in the problem. 
We do not attempt here to calculate the reduction factor 
(which is complicated both technically and conceptually) . 
Instead, we will present the final constraints with the un- 
certainty band corresponding to energy losses as given by 
eq. ^ and reduced by a factor of 5. 

It remains to be checked that, once the PBH becomes 
gravitationally bound, multiple collisions rapidly bring 
the PBH inside the NS. Assuming a radial orbit and de- 
noting the apastron r max , the half-period is 
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The energy loss as a function of r max follows from eqs. ([I]) 
and (|]), 
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where R is the star radius. Dividing the energy loss by 
the time and expressing the energy in terms of r max one 
obtains the differential equation for the evolution of r max 
as a function of time, 
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1 We thank A. Gould for raising this issue. 
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where £ = 



nax /R and 
■kMR 3 ' 2 



r = ~ 6.1 x 10 J s — ^ 

3?tibh V^A^ln A V!0 22 g 

Assuming ^ 3> 1, which is correct for the most part of 
the process, and solving eq. Q the energy loss time is 

£ioss — 2ry^o, 

where the initial value £o can be estimated as £o ~ 
M/ (3tobh hi A), by requiring the initial PBH energy to 
be of the order of Ei om . Assembling all the factors one 
has 
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Thus, PBHs heavier than topbh ^ 3 x 10 17 g end up 

yr. 



inside the NS in a time shorter than 10 10 



B. Capture Rate 

In order to calculate the capture rate, we assume that 
the PBHs follow a Maxwellian distribution in velocities 
with the dispersion v, 
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where tt-bh = Pbh/wbH) Pbh being the density of PBHs 
at the star location. It can be expressed in terms of the 
local DM density pdm as follows, 



fipBH 

Pbh — 7^ Pdm- 



DM 
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Following [16, 22 where relativistic corrections have been 
neglected, the capture rate takes the form 
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where 
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is the capture rate assuming PBHs comprise all of the 
DM, and E\ oss is given by eq.Q. 

Two different regimes are possible. In the case when 
the energy loss is small, Ei oss <C tobhw 2 /3, the exponen- 
tial can be expanded and one gets at the leading order 
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In view of eq. ([2| the capture rate is independent of tobh 
in this regime. In the opposite case E\ oss 3> tobh^ 2 /3 
the exponential in eq. ([9| can be neglected and 
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so that the capture rate decreases with increasing TObh- 
In both cases the capture rate is inversely proportional 
to some power of velocity and is thus maximum for sites 
with high dark matter density pdm where the velocity 
dispersion v is small. 



III. CONSTRAINTS 

As previously mentioned, if a NS captures a PBH, the 
accretion of the NS material onto the PBH rapidly de- 
stroys the star. Therefore, observations of NSs imply 
constraints on the capture rate of PBHs which has to be 
such that the probability of the PBH capture is much 
less than one. In view of eq. ^ these constraints trans- 
late into constraints on the fraction of PBHs in the dark 
matter, ^pbh/^dm- 

Given a NS of the age t-^s, the probability of its sur- 
vival is exp(— tj<isF) with F given by eqs. Q and (|9|. 
Requiring that the survival probability is not small leads 
to the constraint 



^pbh < 1 
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DM 
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Depending on the environment where the NS is located, 
F may vary by many orders of magnitude. The most 
stringent constraints come from sites where Fq is high. 
Among such sites, Globular Clusters (GCs) provide the 
best constraints. 

GCs are compact, nearly spherical collections of stars 
scattered over the Galactic halo. They have ages between 
8 to 13.5 Gyr, and as such are the oldest substructures of 
our Galaxy. GCs are made of population II stars, WDs, 
NSs and black holes. A typical GC has an average radius 
of 30 pc, a core radius of 1 pc and a baryonic mass of 
(a few) x 10 5 M Q [33]. 

The DM content of GCs is a matter of an ongoing 
debate. The distribution of metallicity in GCs is bi- 
modal, indicating two subpopulations formed by differ- 
ent mechanisms [53] . The metal-rich GCs are consid- 
ered to be formed during major gas mergers in proto- 
galaxies [2"5H2"g] . These GCs contain very little DM, if 
any. Instead, the metal-poor GCs are thought to have 
been formed in low-mass dark matter halos at very high- 
rcdshift z ~ 10 — 15 [29H34] , Observations of GCs show 
no evidence of DM halos 35 . This is expected as the ha- 
los should have been tidally stripped due to interactions 
with the Galaxy [36]. The DM content would, however, 
be preserved in the cores of such GCs. In support of this 
picture, it has been found in Refs. [3TJ [35] , using high- 
resolution N-body simulations, that many properties of 
simulated GCs with DM halos are similar to those of ob- 
served GCs. In what follows we will focus on metal-poor 
GCs and assume that they have been formed in DM halos 
and thus possess DM-rich cores. 

In Ref. [37] the DM density close to the GC core 
has been estimated to be of the order prjM ~ 2 x 
10 3 GeVcm~ 3 . Such a DM density was concluded to be 
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FIG. 1. Constraints on the fraction of PBHs in the to- 
tal amount of DM. Shaded regions are excluded. PBHs with 
masses dibh < 10 16 g would have evaporated by now or over- 
produced 7-rays. EROS microlensing survey puts limits on 
the contribution to DM consisting of PBHs in the mass range 
10 26 g < m B H < 3 x 10 34 g from microlensing. In [14] the 
amount of DM adiabatically contracted during star forma- 
tion constrains PBHs in the mass range 10 16 g < mBH < 
3 x 10 22 g. The analysis of [13] constrains the mass range 
10 33 g < tobh < 10 40 g from measurements of the cosmic 
microwave background by FIRAS and WMAP3. The green 
shaded region is excluded by capture of PBHs by NSs in GCs. 
The green hatched region corresponds to constraints that can 
be reached assuming a less conservative choice of parameters 
and a better efficiency of the dynamical friction. 



rather independent of the original halo mass and is in 
agreement with N-body simulations [3TJI3S]- Therefore, 
we adopt this value in our estimates. 

The velocity dispersion is another important parame- 
ter. Since stars are collisionless and therefore behave sim- 
ilarly to DM particles, this parameter can be extracted 
from observations. We adopt the value of v — 7 km s _1 
(the velocity dispersion varies noticeably from cluster to 
cluster; for individual values see Ref. 38 ). Finally, we 
adopt the NS radius i?NS = 10 km, the mass M N s = 
1 M G , and the life time i N s = 10 10 yr [55] . 

The constraints coming from observations of NSs in 
GCs, as well as previously existing constraints are sum- 
marized in Fig. [I] As mentioned earlier, the lack of 
knowledge of the efficiency of the dynamical friction in 
case of NS is reflected in the uncertainty band. As one 
can see, the most conservative constraints exclude the 
PBHs as the unique DM component for masses lower 
than 771bh ~ 5 x 10 23 g, thus pushing up the previously 
existing threshold by more than an order of magnitude. 

In qualitative terms, the shape of the exclusion re- 
gion in Fig. [I] is easy to understand from eqs. ( 10 ) and 
(111. The horizontal part of the curves is due to eq. ( 10 ) 
where the dependence on the PBH mass cancels out (cf. 
eq. (2|)- The inclined part on the right results from 
eq. (11). The transition between the two regimes oc- 
curs at the PBH mass such that E\ oss ~ mQnv 2 /3. The 
left part of the constraints curves comes from the time 



needed for multiple collisions to bring the PBH inside the 
NS. 



IV. CONCLUSIONS 

We have derived constraints on the fraction of PBHs 
in the total amount of DM from observations of NSs in 
globular clusters. These constraints arise from the re- 
quirement that PBHs are captured by NSs with proba- 
bility much less than one, since capture of even a single 
PBH leads to a rapid accretion of the star matter onto 
the PBH and eventual star destruction. 

The most conservative constraints exclude PBH as the 
only DM candidate for masses 10 18 g < tobh < 5 x 10 23 g, 
but accounting for uncertainties can result into the ex- 
clusion of masses 3 x 10 17 g < ttt-bh < 2 x 10 24 g. Taking 
into account the previously existing constraints, only a 
small window of masses around 10 25 g remains where 
PBHs can still constitute all of the DM. Note however 
that this would require a PBH production scenario ca- 
pable of creating a very narrow PBH mass distribution 
function. Those scenarios where the PBH mass function 
is wider than about two orders of magnitude are alto- 
gether excluded. 

The present constraints are complementary to the con- 
straints derived in Ref. 14J. An extra window of masses 
3 x 10 22 g < mBH < 5 x 10 23 g has been constrained at the 
percent level. In Ref. [14] better constraints are achieved 
for masses 10 16 g < ttibh < 8 x 10 19 g, while the present 
paper reaches more competitive constraints for masses 
mBH > 8 x 10 19 g. It is also important to note that dif- 
ferent assumptions are required in these cases: while the 
constraints of Ref. [TJ] are sensitive to the DM distribu- 
tion at the epoch of the GC formation, for the constraints 
derived in this paper the present-epoch DM distribution 
in GCs is relevant. 

We did not present the constraints that come from 
observations of the Galactic center, which is another rel- 
atively close region of high DM density. If the DM den- 
sity in the Galactic center is comparable to that assumed 
above for the the cores of the GCs, no new constraints 
arise from that region [JD]. This is because the capture 
rate depends strongly on the PBH velocity dispersion, cf. 
eq. (10), which is by more than an order of magnitude 
larger in the Galactic center than in the cores of GCs. It 



has been suggested, however, that the DM density in the 
Galactic center may be as high as pdm = 10 6 GeV cm -3 
|41) . If this were confirmed, the constraints from the 
Galactic center would close the remaining allowed PBH 
mass window. 
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